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We model the spectral energy distribution (SED) from the mid-infrared to submil- 
limeter of the ring-like disk of HR4796A, the dustiest A-type star. We consider dust 
made either of coagulated but otherwise unaltered protostellar interstellar grains, or 
grains that are highly-processed in a protostellar /protoplanetary nebula with silicate 
dust annealed and carbon dust oxidized. Both types of dust are successful in repro- 
ducing the observed SED, provided that the grains are highly fluffy, with a vacuum 
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volume fraction of ~ 90%. We find no evidence for the existence of a hot "zodiacal 
dust" component a few AU from the star, which was suggested by previous workers to 
account for the 10 /^m-wavelength emission. 

Subject headings: circumstellar matter — dust, extinction — infrared: stars — plane- 
tary systems: protoplanetary disks — stars: individual (HR 4796A) 



1. Introduction 

HR4796A is a nearby (distance to the Earth d ~ 67 ± 3pc) young main-sequence (MS) star 
(age « 8± 3Myr) of spectral type AO V (effective temperature T eff « 9500 K). The HR4796A dust 
disk has recently aroused considerable interest because (1) it has the largest fractional infrared 
(IR) luminosity relative to the total stellar luminosity (Ltr/L* ~ 5 x 10~ 3 ; L* ~ 21 L@) among 
the ~ 1500 A-type MS stars in the Bright Star Catalogue (Jura 1991); (2) unlike the majority of 
circumstellar disks, the HR4796A disk displays a ring-like structure peaking at ~ 70 AU from the 
central star and abruptly truncated both interior and outside with a width of < 17 AU (Schneider 
et al. 1999); (3) its young-age nature places it at a somewhat transitional stage between massive 
gaseous protostellar disks around young pre-MS T-Tauri and Herbig Ae/Be stars (~ lMyr) and 
much evolved and tenuous debris disks around MS "vega-type" stars (~ 100 Myr). 

The physical and chemical properties of the dust in the HR4796A disk, which play a significant 
role in understanding the creation, dynamical evolution and structural properties of the disk, are 
poorly constrained: (1) Based on the lifetime of dust against loss by Poynting-Robertson drag, 
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Jura et al. (1993) suggested that the HR4796A dust is larger than 10 /mi in radius and the dust is 
a remnant of the protostellar nebula; (2) From an analysis of radiation pressure vs. gravitational 
attraction, Jura et al. (1995) found that the minimum radius of grains orbiting HR4796A is about 
3 /jm and argued that these grains have undergone measurable coalescence; (3) In order to explain 
the > 40 AU inner hole (Jura et al. 1995) and the ss 110 K black-body approximation for the 10- 
100 /jm dust emission, Jura et al. (1998) proposed that the HR4796A grains are largely composed 
of ice particles with a typical radius of ~ 100 /im which originate from a protocometary cloud; (4) 
In order to explain the entire 12.5 /im emission and the 20.8 /tm residual emission unaccounted for 
by their simple model, Koerner et al. (1998) proposed the existence of a tenuous dust component 
at a distance of 3-6 AU from the star which is heated to about 200-300 K, similar to the zodiacal 
dust in our own solar system; (5) From the red reflectance of the disk, Schneider et al. (1999) 
argued that the mean dust size must be larger than several microns and the dust is circumstellar 
debris rather than interstellar in origin; (6) Based on a detailed modelling of the spectral energy 
distribution (SED) of the HR4796A disk, Augereau et al. (1999) proposed a two-component model 
consisting of a cold annulus peaking at 70 AU from the star (made of interstellar dust-type grains 
ranging from 10 /im to a few meters) and a hot population at 9-10 AU from the star (made of 
comet-like grains of radii « 450 /im); (7) Using the flux ratio at 10.8 /im and 18.2 /im and assuming 
"astronomical silicates" , Telesco et al. (2000) inferred a "characteristic" diameter of « 2 — 3 /im 
and argued that these grains are unlikely to be primordial; instead, they are probably products of 
recent collisions of large bodies. 

It is the purpose of this work to constrain the dust properties (size, composition, and morphol- 
ogy) of the HR4796A disk. In general, in a dusty system, dust spatial distribution and dust sizes 
cannot be uniquely determined by the SED alone. Given that the distribution of the dust in the 
HR4796A disk is well constrained by the near-IR imaging of scattered starlight (Schneider et al. 
1999) and, to a less degree, by the mid-IR imaging of dust thermal emission (Jayawardhana et al. 
1998; Koerner et al. 1998; Telesco et al. 2000), an attempt is therefore made in this work to infer 
the HR4796A dust properties by modelling the full spectral energy distribution from the mid-IR 
to the submillimeter wavelengths. 

Lacking a priori knowledge of the composition of the dust in the HR4796A disk, we consider 
two extreme dust types: (1) "cold-coagulation" dust - dust in protoplanetary disks is formed 
through cold aggregation of unaltered interstellar materials - this is the case in the outer parts of 
the disk where interstellar dust originating from the parent molecular cloud (out of which the disk 
forms) is not likely to have significantly changed its composition during passage through the weak 
shock front when the disk is first created, except volatile ice mantles around the refractory cores 
may partly sublimate, but the recondensation of volatiles occurs efficiently behind the shock (see 
Beckwith, Henning, & Nakagawa 2000 for a review); (2) "hot-nebula" dust - dust has undergone 
significant destruction and modification in protoplanetary accretion disks through annealing of 
amorphous silicates and oxidation of carbonaceous dust in the warm inner regions of the accretion 
disk (Gail 2001); turbulent radial mixing (Gail 2001) and/or outflows driven by "X-winds" (Shu, 
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Shang, & Lee 1996) then carries these heavily processed dust into cold regions of the disk where 
it is mixed with freshly accreted material from the parent molecular cloud; as an extreme, we 
consider a dust model in which the carbonaceous dust component has been fully destroyed by 
oxidation so that only crystalline silicate dust remains. In reality, the dust in protoplanetary disks 
would undoubtedly be intermediate between heavily processed dust and aggregates of unaltered 
interstellar dust. Since planetesimals and cometesimals are formed by coagulation of such dust 
aggregates, the dust generated by collisions of planetesimals and cometesimals should resemble the 
original dust from which they are built up. 

We first discuss in §2.1 the general constraints on the HR4796A dust and disk properties. 
We then model the observed SED in §2.2 and §2.3 in terms of the "cold-aggregation" model of 
unaltered interstellar dust and the "hot-nebula" model with all amorphous silicate dust annealed 
and all carbonaceous dust oxidized, respectively. We find that, for both models, (1) highly porous 
grains with a porosity P « 0.90 provide excellent fits to the entire SED; (2) there appears to be no 
need for a hot "zodiacal cloud" dust component (§3.3) which was previously suggested by Koerner 
et al. (1998) and Augereau et al. (1999). Therefore, the presence of either type of particles in the 
disk is possible. We show in §3.2 that, although good fits to the SED can be achieved by models 
with a single power-law dust spatial distribution, they are apparently in conflict with the imaging 
observations of scattered starlight and dust thermal emission. In §3.5 we calculate the dust IR 
intensities integrated over the SIRTF/IRAC and MIPS bands predicted for our best-fitting models. 
The major conclusions are summarized in §4. 

2. Modelling the Dust IR Emission 

A grain in the optically thin dust disk of HR 4796A absorbs stellar ultraviolet / optical photons 
and then re-radiates the energy in the IR. To calculate the disk's emission spectrum, knowledge of 
its morphology, composition, and size is required. 

2.1. General Considerations on Dust and Disk Properties 

A fluffy structure is expected for dust in protoplanetary disks as a result of its coagulational 
growth process. We characterize its fluffiness by porosity P, the fractional volume of vacuum. We 
assume all grains are spherical in shape. 

In §A we estimate that the "cold-coagulation" model leads to porous dust consisting of 
amorphous silicate and carbonaceous materials (and H^O-dominated ices in regions colder than 
~ 110 — 120 K) with a mixing ratio of m car b/w- s ii ~ 0.7 [and m- lce /{m s \\ + m car b) ~ 0.8 for cold 
regions]. For the "hot-nebula" dust model, we assume that the dust is exclusively composed of 
crystalline silicate dust (and ices in cold regions). 
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We assume a power law dust size distribution dn(a)/da oc a~ a which is characterized by a 
lower-cutoff a m i n , upper-cutoff a max and power-law index a (where a is the spherical radius). We 
take a m i n = 1 fim, the smallest value among previous estimations of a m i n (see §1);* and a max = 1 cm 
(this is not a critical parameter since grains larger than ~ 100 /jm are like blackbodies and their IR 
emission spectra are size-insensitive). 

The dust spatial distribution is well constrained by the near-IR imaging of scattered starlight 
(Schneider et al. 1999) to be a sharply-truncated ring-like structure peaking at a radial distance 
of ~ 70 AU from the star with a characteristic width of < 17 AU. Therefore, we adopt a Gaussian 
function for the dust spatial density distribution: 2 dn{r)/dr oc exp[— 41n2{(r — r p )/A} 2 ] which 
is characterized by the radial position r p where dn(r)/dr peaks and the full width half maximum 
(FWHM) A. The dust surface density distribution obtained by integrating dn/dr over the perpen- 
dicular (to the disk plane) path lengths can be expressed as a(r) = <7 p exp[— 41n2{(r — r p )/A} 2 ] 
where a p is the mid-plane (z = 0) surface density at r = r p . This distribution function, with 
r p = 70 AU and A = 12 AU, was successful in modelling the scattered-light images (Kenyon et al. 
1999). A Gaussian- type dust distribution was also adopted by Klahr & Lin (2000) to model the 
dynamics of the dust ring around HR4796A. We fix r p to be at r p = 70 AU and take A = 15 AU, 
the mean value of the determinations of Schneider et al. (1999; A < 17 AU) and Kenyon et al. 
(1999; A = 12 AU). The inner boundary r Ui is set to be where grains are heated to > 1500 K, 
and hence, r Ui is a function of dust size. For micron-sized grains, the inner boundary is roughly 
r- m = 0.15 AU. The outer boundary is taken to be r out = 250 AU which is expected from the disk 
truncation caused by the tidal effects of HR 4796B , a companion star of HR 4796A ( Jayawardhana 
et al. 1998). 

Therefore, we are only left with two free parameters: (1) the dust porosity P, and (2) the dust 
size distribution power index a. 

Dielectric functions are taken from (1) Draine & Lee (1984) for amorphous silicate dust; (2) Li 
& Draine (2001a) for crystalline silicate dust; (3) Li & Greenberg (1997) for carbonaceous dust; (4) 
Li & Greenberg (1998) for H^O-dominated ice. The Bruggman effective medium theory (Bohren 
k, Huffman 1983) is used to calculate the mean dielectric functions for the fluffy heterogeneous 
dust aggregates. Absorption cross sections are obtained using Mie theory. Approximating the 
HR4796A radiation field by the Kurucz model atmosphere spectrum for AO V stars (Kurucz 1979), 
dust equilibrium temperatures are then derived by balancing absorption and emission. For a given 



X A fluffy grain of a = 1 fim with a porosity of P = 0.90 consists of ~ 100 constituent individual (interstellar) 
particles which have a typical size of a ~ 0.1 fxra (see Li & Greenberg 1997). Models with a smaller a m i n (= 0.1 fim) 
and a larger a m i n (= 10 /jm) will be discussed in §3.2 

2 Following Kenyon et al. (1999), we assume an exponential, radial-independent, vertical distribution dn/dr oc 
exp[— (z/\/2~ff) ] where the vertical scale height H is ~ 0.5 AU. Alternatively, H can be determined from the 
vertical hydrostatic equilibrium assumption (H oc r 3//2 ; see Appendix B in Li & Lunine 2003). But the knowledge 
of the dust vertical distribution is not required and it does not affect our results since in modelling the SED of the 
HR4796A disk, what actually involves is the dust surface density distribution er(r). 



5 



dust size distribution and a given disk structure (dust spatial density distribution), the emergent 
IR emission spectrum can be obtained by integrating over the dust size range, and over the entire 
disk. The calculated IR spectrum is then compared with the available photometric data for the 
HR4796A disk compiled by Augereau et al. (1999). 

2.2. "Cold-Coagulation" Dust Model 

We first model the HR4796A SED in terms of the dust generated from collisions of plan- 
etesimals and cometesimals formed in the disk (see §3.4) through cold-coagulation of unaltered 
interstellar grains from its parent molecular cloud. Using the dust composition discussed in §A, 
we calculate the model IR spectra for a wide range of dust porosities P and a wide range of dust 
size distribution power-indices a to search for good fits. We note that P refers to the porosity of 
refractory dust - P will be reduced to P' for dust with ices (see §B). 

We illustrate in Figure 1 the best fits obtained for dust with a porosity of P = 0.95, 0.90, 0.80, 0.60. 
Model parameters and results are tabulated in Table 1. It is seen that the best fit is provided by 
dust with P ~ 0.90 (with a total mass of f» 0.67 m ffi ; model no. 2) - dust more porous than this 
is somewhat too hot so that its emission is deficient in the submillimeter wavelength range (see 
Figure la); on the other hand, dust more compact than this is a bit too cold so that it produces 
too much emission in this wavelength range (see Figure lb). 

2.3. "Hot-Nebula" Dust Model 

Now we consider another extreme case: dust in the inner disk regions has been so heavily 
processed that all amorphous silicate dust has been annealed (T > 800 K) and all carbonaceous dust 
has been oxidized (T > 1100 K) by reacting with OH. The annealed silicate dust is then transported 
to (or the silicate vapour recondenses in crystalline form in) the cold outer regions of the disk where 
these grains grow into fluffy aggregates and are ultimately built into planetesimals/cometesimals. 
Apparently, only a fraction of the dust in the disk is accreted to the inner warm regions and, hence, 
only a fraction of the silicate dust is crystalline. But since we are considering an "extreme" case, in 
this section we model the "hot-nebula" dust as porous aggregates of pure crystalline silicate dust 
(plus ices in cold regions where dust reaches a temperature of ~ 110 — 120 K). We do not know how 
much ice would recondense on the silicate core seeds in the fluffy aggregate, but the assumption of 
full condensation of all condensible volatile elements (C, O, N) as ices (see §A) seems to be at the 
highest end. As a complement to §2.2 where full condensation is assumed, we take Vi ce /v s \\ = 1, a 
plausible value for dense molecular clouds. 

Again, model IR spectra are calculated for a wide range of dust porosities P and a wide range 
of dust power-indices a. Except the sharp features at 11.3 /um and 23 /xm, the "hot-nebula" model 
spectra are very similar to those of the "cold-coagulation" model. For illustration, we show in Figure 
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2 the best-fit model spectra provided by dust with P = 0.95, 0.90, 0.80, 0.60. Model parameters and 
results are also tabulated in Table 1. Similar to the cold-coagulation model, the best fit is given by 
dust with P ~ 0.90 (with a total mass of ps 1.25 m®; model no. 6). 



Table 1: Models for the HR4796A dust disk IR emission^ 
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1 


cold-coag. 


Gaussian 


0.95 


0.87 


15 




2.8 


2.25 


47.8 


9.46 


2.67 


2.97 


0.041 


2.85 




2 


cold-coag. 


Gaussian 


0.90 


0.73 


15 




2.9 


2.11 


28.7 


4.34 


4.03 


4.90 


0.042 


1.81 


preferred 


3 


cold-coag. 


Gaussian 


0.80 


0.46 


15 




3.0 


2.00 


18.4 


2.00 


5.73 


7.55 


0.043 


6.08 




4 


cold-coag. 


Gaussian 


0.60 





15 




3.2 


1.83 


9.26 


0.44 


4.34 


13.1 


0.038 


7.60 




5 


hot-ncbu. 


Gaussian 


0.95 


0.90 


15 




2.6 


2.66 


155 


45.5 


6.02 


0.85 


0.040 


5.29 




6 


hot-ncbu. 


Gaussian 


0.90 


0.80 


15 




2.8 


2.25 


47.8 


9.46 


7.48 


2.55 


0.036 


4.48 


preferred 


7 


hot-ncbu. 


Gaussian 


0.80 


0.60 


15 




2.9 


2.11 


28.7 


4.34 


11.6 


4.30 


0.037 


5.24 




8 


hot-ncbu. 


Gaussian 


0.60 


0.20 


15 




3.1 


1.91 


12.6 


0.93 


9.25 


8.04 


0.031 


9.94 




9 


cold-coag. 


power-law 


0.90 


0.73 




-2.3 


3.2 


1.83 


9.26 


0.44 


8.43 


0.42 


0.0011 


2.32 




10 


cold-coag. 


power-law 


0.90 


0.73 




-2.5 


3.0 


2.00 


18.4 


2.00 


5.81 


0.98 


0.0054 


3.28 


r S [40, 13 


11 


hot-ncbu. 


power-law 


0.90 


0.80 




-1.1 


2.9 


2.11 


28.7 


4.34 


25.5 


0.26 


0.0022 


9.71 




12 


hot-ncbu. 


power-law 


0.90 


0.80 




5.2 


2.6 


2.66 


155 


45.5 


6.29 


0.056 


0.0026 


2.02 


r e [40, 13 


13 


cold-coag. 


Tclcsco 


0.90 


0.73 


21 




2.9 


2.11 


28.7 


4.34 


5.41 


2.34 


0.020 


3.60 




14 


hot-ncbu. 


Tclcsco 


0.90 


0.80 


21 




2.8 


2.25 


47.8 


9.46 


10.0 


1.22 


0.017 


4.89 




15 






0.90 


0.73 


15 






0.22 


0.81 


0.15 


5.49 


193 


0.045 


3.44 


a min = 0. 


16 


cold-coag. 


Gaussian 


0.90 


0.73 


15 




3.1 


19.1 


1048 


117 


2.98 


0.13 


0.042 


5.88 


a min = 10 


Unless otherwise stated, all models 
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ttmin 


= 1 fim, a max = 1 cm, r p = 


70 AU, and 


r € [0.15, 250 AU]. 








Dust spatial distribution - "Gaussian": o(r 
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p exp[ 


-41n2{(r- 


r- P )/A} 2 ]; 


"power-law" 


: a(r) 


= <T P (r/r p y 


"Tele- 






sco": a 


(r) = <7 P exp 


[-(A/AU) In 2 (r/r p )] for r 


e [0.15, 70 AU]; a 


(r) = ffp/{l + [(r-r P 


)/A] 2 } for r £ [70, 105 AU]; 







<r(r) = op [-0.01 (r/AU) + 1.32] for r e [105, 130 AU]. 

c (a 7 } = J' 0max da a 1 dn(a) /da/ J^ max da dn{a) /da. The dust surface mass density can be written as a(r; m) = (m) a(r) 
with the mean dust mass (m) = (47r/3)(a 3 )(p) where (p) w 2.5(1 - P)gcm~ 3 , (p'} w 1.7(1 - P')gcm -3 for the 
"cold-coagulation" model; and (p) w 3.5(1 — P)gcm~ 3 , (p 1 ) w 2.4(1 — P')gcm~ 3 for the "hot-nebula" model (see 
§B). 

d Vertical optical depth at A = 0.55 pm and r — r p . The reason why r p is much smaller for models no. 9-11 is that 
the bulk of the dust is piled at the edge of the disk which is unphysical (see text [§3.2]). 



3. Discussion 

It is seen in §2 that, using a simple power-law dust size distribution {dn/da oc a~ a with 
a ~ 2.8 — 2.9, (a) ~ 2 fj,m) and a Gaussian-type dust spatial distribution with a peak at r p = 70 AU 
and a FWHM A = 15 AU as inferred from the near-IR imaging observation of scattered starlight 
(Schneider et al. 1999; Kenyon et al. 1999), both the "cold-coagulation" dust model (no. 2) and 
the "hot-nebula" dust model (no. 6) are successful in reproducing the observed SED, provided that 
the dust in the HR4796A disk is very fluffy, with a porosity P w 0.90. We note that, since we use 
this dust spatial distribution, it is expected that our models are also able to reproduce the imaging 
observations. 

Both the "cold-coagulation" model and the "hot-nebula" model predict a vertical optical depth 
at visible wavelengths of t v ~ 0.04 at r = r p (see Table 1). This justifies the optically-thin 
treatment employed in the entire paper. This is because that in these models, with a power-law 
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index for the size distribution of a G [2.6,3.2] (see Table 1), most of the mass of the particles is 
in the largest grains. However, the opacity per gram dust (ft a b s = [3/4p] Qabs/ a where Q a bs is 
the absorption efficiency) for these macroscopic grains (Q a bs ~ 1 at visible wavelengths) is much 
smaller than that of micron-sized dust: K a b s oc 1/a. We note that although the maximum dust 
size a max is not well constrained, the predicted IR emission spectrum and the optical depth are not 
very sensitive to the precise value of a max since grains larger than ~ 100 fim have Q a b s ~ 1 and 
emit like black-bodies. 



3.1. Robustness 

The extent to which dust in protostellar disks has been processed prior to incorporation into 
planetesimals and comets is not known. But it seems plausible that a fraction of precometary 
materials has been processed in protostellar nebulae, based on the observationally presence of 
crystalline silicate dust in comets (see Wooden 2002 for a summary of the evidence). While at 
most ~ 5% of the silicate material is crystalline in the interstellar medium (Li Sz Draine 2001a), 
in some comets the fraction of crystalline silicate dust is as high as ~ 30%. Accretion heating in 
protoplanetary disks is a probable source for the formation site of the crystalline silicate - although 
26 Al heating of comets and consequent release of heat in the amorphous-to-crystalline ice transition 
cannot be ruled out. The warm inner zone of the accretion disk is particularly appealing as a site 
for the annealing of amorphous silicate. As a result of diffusional mixing induced by turbulence in 
the optically thick disk (Gail 2001), or outflows driven by reconnecting magnetic field lines (the 
"X-winds" model of Shu, Shang, & Lee [1996]), a fraction of the annealed silicate dust arrives at 
the cooler outer part of the disk to be incorporated in icy "cometary" bodies. 

Therefore, we expect that the dust generated by collisions among cometary bodies (see §3.4), 
for either origin of the crystalline silicate material, should be intermediate between the "cold- 
coagulation" dust and the "hot-nebula" dust. The mixing ratio of these two dust types is unclear. 
In principle, the 10 /im silicate emission feature should allow us to infer the thermal history of 
silicate dust through the absence or presence (and the mass fraction) of crystalline silicate dust. 
However, the only available 8-13 /im spectrum shows very weak thermal emission in the silicate 
feature; its quality is inadequate to discern the detailed spectral features expected from crystalline 
silicate dust (Sitko, Lynch, & Russell 2000). But in any case, the fact that the two extreme models 
are able to provide close fits to the observed SED, any linear combinations of these two dust types 
(i.e., any degree of processing/modification) are expected to be able to reproduce the observed SED 
as well. Therefore, our models are robust. 
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3.2. Uniqueness 

In modelling SED, dust sizes and dust spatial distribution are coupled. Since the dust spatial 
distribution we adopted for the HR 4796A disk is well constrained by near-IR images of scattered 
light (Schneider et al. 1999), the dust properties are more or less uniquely determined: the dust 
must be very fluffy (with P ~ 0.9); the dust power-law size distribution has an index a pa 2.8 — 2.9 
with a mean dust size of ~ 2 fj,m (see Table 1); the dust composition should be somewhere between 
the "cold-coagulation" dust and the "hot-nebula" dust although its detailed composition (e.g. the 
crystallinity of the silicate dust) is yet unknown (see §3.1). 

The dust spatial distribution in the HR4796A disk was often modelled as a single power-law 
(Jayawardhana et al. 1998; Koerner et al. 1998; Wyatt et al. 1999) or two power-law components 
(Augereau et al. 1999). For the purpose of comparison, we have also modelled the HR4796A SED 
assuming a single power-law dust distribution dn/dr oc r~@ in terms of the P = 0.90 (P' pa 0.73) 
"cold-coagulation" dust model and the P = 0.90 (P' pa 0.80) "hot-nebula" dust model. 3 We 
consider two disk extents: (1) r m = 0.15 < r < r out = 250 AU (see §2.1); and (2) n n = 40 < 
r < r out = 130 AU. The latter takes into account the existence of an inner hole with a radius of 
~ 40 — 60 AU which was first noticed by Jura et al. (1995) from an analysis of the IRAS (Infrared 
Astronomical Satellite) and ground-based photometry. The inner hole was later confirmed by mid- 
IR imaging of dust thermal emission carried out by three independent groups (Jayawardhana et 
al. 1998; Koerner et al. 1998; Telesco et al. 2000). These mid-IR imaging observations also implied 
an outer radius of ~ 110 — 130 AU for the HR4796A disk. As can be seen in Figure 3, models 
with a single power-law dust spatial distribution are also successful in reproducing the observed 
SED. However, except the 40 < r < 130 AU "hot-nebula" model (no. 12), all models (no. 9 - no. 11) 
require an increasing accumulation of dust toward the outer edge (i.e. (3 < 0). This is neither 
physical nor consistent with the imaging observations of scattered light (Schneider et al. 1999) and 
dust thermal emission (Jayawardhana et al. 1998; Koerner et al. 1998; Telesco et al. 2000). 

It is expected that models with two power-law dust spatial distributions are also able to 
reproduce the observed SED. But the constraints placed by the near-IR and mid-IR imaging ob- 
servations would imply that the two-segment power-law distribution should not deviate much from 
the Gaussian function. 

Assuming spherical silicate dust of a diameter of 2.5 yum (estimated from the 10.8 /jm and 
18.2 fim emission), Telesco et al. (2000) derived the dust spatial distribution from the 18.2 /im 
brightness distribution (see their Figure 5). We approximate their result as: for r < r p = 70 AU, 
a pa dp exp [- (A/AU) In 2 (r/r p )] (A pa 21 AU); for 70 AU < r < 105 AU, a pa a p / {l + [(r - r p ) /A] 2 
for 105 AU < r < 130 AU, a Pa a p [-0.01 (r/AU) + 1.32]. In Figure 4 we show that the best-fitting 
"cold-coagulation" (no. 13) and "hot-nebula" models (no. 14) with this dust distribution also pro- 
vide close fits to the observed SED. 



3 As in §2.1, the dust surface density distribution is written as a(r) = a p (r/rp) 
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The overall shape of the Telesco et al. (2000) distribution resembles that of the Gaussian 
distribution, except that the former is relatively broader and has a flatter wing at r > 70 AU. 
Since the Gaussian distribution was derived from scattered light images which reflect dust of all 
sizes while the Telesco et al. (2000) distribution was derived only from the 10.8 /mi and 18.2 /mi 
mid-IR-emitting dust, we prefer the Gaussian distribution although the real dust distribution may 
not be an exact Gaussian function. 

To illustrate the disk regional contribution, we show in Figure 5 the IR emission from dust 
in the regions of (r p — A) < r < r p , r p < r < (r p + A), and (r p — A) < r < (r p + A) as well 
as 0.15 AU < r < 250 AU calculated from the best-fitting "cold-coagulation" dust model (no. 2; 
P = 0.90, P' « 0.73, a « 2.9) with a Gaussian spatial distribution (r p = 70 AU, A = 15 AU). 
It is seen that the IR emission is exclusively contributed by dust confined in the ring region of 
(r p -A) <r < (r p + A). 

So far, all models assume a m j n = 1 /mi. We now consider models with smaller or larger ftmin 
values. In Figure 6a we show the best-fit spectra calculated from models with o m i n = 0.1 /im (model 
no. 15) and a m - m = 10 /mi (model no. 16). It is seen that (1) models with a m j n = 0.1 ^m 4 emit a 
bit too much at A > 100 p; decreasing the dust porosity improves the fit at A > 100 /im but the 
fit to the A < 10 /im part deteriorates; (2) models with a m \ n = 10 /im emit too little at A < 10 /im; 
increasing the dust porosity does not solve this problem. 

To illustrate the IR emission contributed by different dust sizes, we show in Figure 6b the 
emission spectra calculated from the best-fitting "cold-coagulation" model (no. 2; P = 0.90, P' ps 
0.73, a ~ 2.9) with 1 /im < a < 10 /mi (subject to radiative expulsion; see §3.4), 1 /mi < a < 100 /mi 
(subject to Poynting- Robertson inward spiralling drag; see §3.4), and 1 /im < a < 1cm. It is seen 
that the A < 20 /im mid-IR emission is dominantly produced by grains smaller than 100 /xm in 
radius. Since these grains will be removed from the disk during its lifetime either through radiative 
expulsion or through Poynting-Robertson drag, there must exist a source of replenishment (see 
§3.4). 

Finally, we show in Figure 7a the IR emission per unit dust mass received at the Earth 
produced by the P = 0.90 (P f « 0.73) "cold-coagulation" dust at r = 70 AU with a single size of 
a = 1 /xm, 10 /im, 100 /im, 1000 /im, 1 cm (at r = 70 AU, all grains are colder than 110 K and have a 
porosity of P' = 0.73 [their constituent particles are coated by a layer of ice mantles] except the 
a = 1/mi dust which attains an equilibrium temperature of ~ 154 K and thus has no ice mantles 
and has a porosity of P = 0.90). In Figure 7b we show the IR spectra emitted by a gram of the 
P = 0.90 (P' » 0.73) "cold-coagulation" dust at r = 55 AU, 70 AU, 85 AU, 130 AU, 250 AU with a 
distribution of sizes (dn/da cx a~ a , a = 2.9, a € [1 /mi, 1cm]). 



4 This requires that the constituent individual particles should be smaller than the typical interstellar dust (a ~ 
0.1 /im). 
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3.3. Is There A "Zodiacal Dust" Component? 

In modelling the HR4796A SED, Koerner et al. (1998) and Augereau et al. (1999) argued that, 
in addition to the dust responsible for the emission at A > 20 /im, a hot "zodiacal dust" component 
which is confined within ~ 3 — 6 AU or ~ 9 — 10 AU of the star is required to account for the entire 
emission at A < 12 /im and < 5 — 10% of the emission at A ~ 20 /im. However, it can be seen clearly 
from Figure 1, Figure 2, and even the unphysical single power-law spatial distribution model of 
Figure 3 that, our models closely reproduce the observed SED over the entire wavelength range. 
There is thus no need to invoke a hot "zodiacal" dust component; indeed we may now ask how 
much such dust can be tolerated based on our models. 

An upper limit on the "zodiacal" dust component can be obtained such that this abundance of 
hot dust, after added to the IR emission spectra of our best-fitting models (e.g., no. 2), would not 
exceed the observed emission at A < 12 /tin. Following Augereau et al. (1999), we assume that the 
"zodiacal" dust consists of porous (P « 0.97) crystalline silicate grains of a = 450 /tni. As shown 
in Figure 8, if the observed emission at A < 12 /im is entirely produced by the "zodiacal" dust, a 
total amount of rrid ~ 9.5 (1.8) x 10 -5 is required if the "zodiacal" dust lies at r = 9 (4.5) AU. 
It is seen in Figure 8 that our best-fitting model (no. 2) allows at most 10% of this mass to be in 
the "zodiacal" dust component. 

The reason why Koerner et al. (1999) required a population of "zodiacal dust" can be ascribed 
to the fact that they have not considered a distribution of dust sizes; instead, they simply approx- 
imated the dust emissivity e\ as a one-parameter function: ex = 1 for A < ao and e\ = ao/A for 
A > ao where ao is a parameter characteristics of dust size. It would not be surprising if their 
model lacks hot dust emitting at A < 12 /im since their approximation is essentially for single-sized 
dust models. 

The reason why Augereau et al. (1999) needed a population of "zodiacal dust" may lie in the 
fact that they fixed the power-index of their dust size distribution to be a = 3.5. Using a steeper 
size distribution (in comparison with a = 2.8 — 2.9 of our P = 0.90 best-fitting models [no. 2 and 
no. 6; see Figures 1,2], Augereau et al. (1999) had to raise a m - m to > 10 /im in order to reproduce 
the long wavelength emission. This led to a paucity of small grains which are hot enough to emit 
at A < 12 /im (this can also be seen in Figure 6a [dotted line]). 

Also arguing against attributing the A < 12 /im emission to a hot "zodiacal dust" component 
are the recent 10.8 //m and 18.2 /on mid-IR images of the HR4796A disk that show that the disk's 
size at 10 /im is comparable to its size at 18/im (Telesco et al. 2000). This implies that the 18/xm- 
emitting dust may also emit some, or all, of the 10 /im radiation (Telesco et al. 2000). This is 
in sharp contrast to the "zodiacal dust" assumption which suggested that the "zodiacal dust" 
component contributes to only < 5 — 10% of the 20 /im emission (Koerner et al. 1998; Augereau et 
al. 1999). 

The tidal effects by sweeping inner planets together with the Poynting-Robertson drag may 
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be responsible for the absence of a warm "zodiacal dust" population as well as the presence of an 
inner disk hole (r < 50 AU). 

3.4. Radiation Pressure and Poynting-Robertson Drag 

In addition to the gravitational attraction from the central star, grains in the HR4796A disk 
are subject to (1) radiative repulsion due to the momentum carried by stellar photons and (2) 
Poynting-Robertson drag which takes both energy and momentum from their orbits and causes 
them to spiral toward the gravitational force center (Burns, Lamy, & Soter 1979; Backman & 
Paresce 1993). 

We have calculated /3rp - the ratio of radiative pressure (RP) force to gravitational force for 
the best-fitting "cold-coagulation" dust (no. 2; P = 0.90, P' ps 0.73) and "hot-nebula" dust (no. 6; 
P = 0.90, P' ps 0.80). As shown in Figure 9a, for grains smaller than ~ 10 — 20 /um, the radiation 
pressure overcomes the gravitational attraction (i.e. /3rp > 1) and, therefore, these grains will be 
blown out from the HR4796A disk. The dust removal rate due to the radiation pressure expulsion 
can be estimated by integrating (47r/3) pa 3 /rRp(r) over the size range 1 /xm < a < 10/xm and 
over the entire disk. Assuming the RP timescale is comparable to the local dynamical timescale 
[trp(t) ps 370 (r/70AU) 3 / 2 yr], we estimate the RP dust mass loss rate to be ps 8.19 x 10~ 7 m©/yr 
for the best-fitting "cold-coagulation" model (no. 2) and ps 8.28 x 10~ 7 m ffi /yr for the best-fitting 
"hot-nebula" model (no. 6). 

We have also calculated the Poynting-Robertson (PR) drag timescales (tpr). In Figure 9b 
we show that, for grains smaller than ~ 100 — 200 jim at a radial distance of r = 70 AU from the 
central star, their lifetimes (at r = 70 AU) due to the PR drag effect are shorter than the HR4796A 
age and hence, for grains in the range of 10 — 20 [xm < a < 100 — 200 fxm, although stable against 
radiation pressure ejection, the Poynting-Robertson drag could remove them from the disk. Since 
the best-fitting models have < 30 — 40% of the total surface areas in grains smaller than 100 /im, 
these grains must be efficiently replenished by cascade collisions of planetesimals and larger grains. 
By integrating the PR dust removal rate (47r/3) /oa 3 /rpR(a, r) over the whole size range and over the 
entire disk, we estimate the PR dust mass loss rate to be ps 9.22 x 10~ 9 m®/yr for the best-fitting 
"cold-coagulation" model (no. 2) and ps 7.96 x 10~ 9 m©/yr for the best-fitting "hot-nebula" model 
(no. 6). 

Therefore, over the life span of HR4796A, roughly 6.7 m® of dust is lost by radiation pressure 
and Poynting-Robertson drag, which is about 3 times smaller than the estimation of Augereau 
et al. (1999; ps 20 m©). This is mainly because Augereau et al. (1999) attributed the dominant 
dust removal process to grain-grain collisions. We consider this unlikely since grain-grain collisions 
would not remove the dust from the disk, instead, they just re-distribute the dust over different 
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size bins through fragmentation or sticking. 5 

It is well established (see e.g. Spangler et al. 2001, Zuckerman & Becklin 1993) that the 
amount of IR excess around main sequence stars decreases with the star's age. There are at least 
two possible explanations for this result: (1) the particles are lost by the Poynting-Robertson drag 
effect; (2) the particles have coalesced into larger bodies which are undetectable as the opacity per 
gram decreases (K a bs oc 1/a, see §3). For the HR4796A disk, the latter process seems to be the 
dominant one as reflected by its collisionally-replenished secondary nature - the amount of dust 
incorporated into larger bodies which act as a source of replenishment must be much larger than 
that of the dust lost through Poynting-Robertson drag. 



3.5. Predictions for SIRTF 

The Space Infrared Telescope Facility (SIRTF) will be capable of sensitive imaging using the 
Infrared Array Camera (IRAC) at 3.6, 4.5, 5.8, and 8.0 /im, and using the Multiband Imaging 
Photometer (MIPS) at 24, 70, and 160 fim. In Table 2 we show the band-averaged intensities for 
our preferred dust models. 

SIRTF will also be able to perform low-resolution 5-40 /im and high-resolution 10-37 /im spec- 
troscopic observations using the Infrared Spectrograph (IRS) instrument. IR spectroscopy and 
imaging will provide powerful constraints on the HR4796A dust spatial distribution and its chem- 
ical composition. The "cold-coagulation" dust models all predict an almost featureless SED (the 
9.7 fim amorphous silicate feature is very broad and smooth). 6 In contrast, the "hot-nebula" dust 
models predict two sharp crystalline silicate features at 11.3 and 23 /im. In order to produce no- 
ticeable 11.3 and 23 /im features, at least 20% of the silicate dust must be in the crystalline form. 
SIRTF spectroscopy will allow us to infer the degree of processing which the HR4796A dust has 
experienced. 



5 Since the grain-grain collision timescale t co ii ~ {2n (a 2 ) a (r)Q(r)} 1 w 670yr at r = 70 AU [fl(r) = (GM*/r 3 ) 1/2 
is the Keplerian frequency; M* (ss 2.5 m©) is the stellar mass; G is the Gravitation constant] is comparable to the 
RP timescale trp « 370 yr (at r = 70 AU), a stable population of a < 10 /im dust can also be obtained through a 
cascade production during collisions between larger grains. 

6 The 3.3, 6.2, 7.7, 8.6 and 11.3 /xm vibrational bands characteristic for polycyclic aromatic hydrocarbon (PAH) 
molecules are not seen in our model spectra although their presence in the HR 4796A disk is expected in the context 
of the "cold-coagulation" model since PAHs are a significant constituent of interstellar dust (see e.g. Li & Draine 
2001b). This is mainly because PAHs have condensed in the icy mantles of the porous aggregates at r > 70 AU. In 
the inner regions of the disk, PAHs, sublimated from the icy mantles, would be rapidly photodissociated. Therefore, 
PAHs are not considered in this paper. The presence of PAHs in dusty disks will be discussed in a subsequent paper. 
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Table 2: Dust IR emission (Jy) averaged over SIRTF bands predicted for our preferred models. 
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IRAC 




IRAC 




IRAC 
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3.6 fim 
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8.0 fim 
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70 /im 
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1.18xl(T 8 
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-7 
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-5 


7.39x10" 


■3 
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7.26 


2.94 


6 


3.44xl0~ 9 


3.91x10" 


-7 


2.18x10" 


-5 


1.67x10" 


-3 


2.85 


6.61 


2.58 


13 


1.31xl(T 8 


9.51x10" 


-7 


4.62x10" 


-5 


6.62x10" 


-3 


2.56 


7.81 


3.41 


14 


3.68xl(T 9 


3.74x10" 


-7 


1.98x10" 


-5 


1.51x10" 


3 


2.79 


7.12 


2.99 



4. Conclusion 

We have modelled the mid-infrared to submillimeter spectral energy distribution of the ring- 
like dust disk around the dustiest ~ 8Myr-old A- type star HR4796A. We start with two extreme 
dust types: one formed through cold-coagulation of unaltered protostellar interstellar grains, and 
one formed through aggregation of grains highly-processed in the protostellar nebula with silicate 
dust annealed and carbon dust oxidized. We adopt a Gaussian-type dust spatial distribution with a 
peak at 70 AU from the central star and a FWHM of 15 AU as inferred from images of scattered light 
and dust thermal emission. We take a simple power-law for the dust size distribution {dn/da oc a~ a 
with q s=s 2.8 — 2.9) in the range of 1 fim < a < 1 cm. Our principal results are: 

1. It is shown that both types of dust are successful in reproducing the observed SED, provided 
that the dust generated by collisions of planetesimals and cometesimals is highly fluffy, with a 
vacuum volume fraction of ~ 90% (§2). The fact that the dust in the HR4796A disk must be 
somewhere intermediate between these two types implies that our models are robust. Future 
high resolution spectroscopy of the HR 4796A disk would allow us to infer the degree to which 
the dust has been processed (§3.1). 

2. Although models with a single power-law spatial distribution are also able to reproduce the 
observed SED, the derived dust distribution is both unphysical and inconsistent with the 
imaging observations of scattered light and dust mid-IR thermal emission (§3.2). 

3. Our models show no evidence for the existence of a hot "zodiacal dust" component suggested 
to lie in a radial distance of a few AU from the star by previous workers to account for the 
emission at A < 12 fim (§3.3). Our upper limit on the total mass of this component is at least 
10 times smaller than what would be obtained if the entire A < 12 fim emission is attributed 
to this hot dust. 

4. Grains smaller than ~ 10 — 20 fim will be radiatively expelled from the disk; grains at r = 
70 AU in the range of ~ 10 — 20 fim < a < 100 — 200 fim will also be removed from the disk 
due to the Poynting-Robertson inward spiralling drag (at a closer distance from the star, 
larger grains will be removed). Collisions of planetesimals/cometesimals must continuously 
replenish the dust in the disk at a rate of « 8.3 x 10~ 7 m e /yr (§3.4). 
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A. Possible Dust Composition for Cold Protoplanetary Disk 

The coagulation of interstellar grains that results in fluffy and inhomogeneous aggregates oc- 
curs in cold, dense molecular clouds and protostellar and protoplanetary dust disks. It plays an 
important role in the formation of planetary systems (Weidenschilling & Cuzzi 1993). In this 
section we approximately derive the proportional composition of the dust in circumstellar disks 
around (pre-)MS stars from the abundances of the condensible elements (C, N, O, Si, Fe, and 
Mg), 7 assuming protostellar activities impose little modification on protostellar grain compositions 
(see Beckwith, Henning, & Nakagawa 2000). 

Let [X/H]q be the cosmic abundance of X relative to H (we assume the cosmic elemental abun- 
dances are those of the solar values: [C/H] Q ss 391 parts per million (ppm), [N/H] Q rj 85.2 ppm, 
[O/H] » 545 ppm, [Mg/H] » 34.5ppm, [Fe/H] « 34.4ppm, and [Si/H] » 28.1ppm [Sofia & 
Meyer 2001]); [X/H] gas be the amount of X in gas phase ([C/H] gas « 140ppm, [N/H] gas « 61 ppm, 
[0/H] gas Ri 310 ppm; Fe, Mg and Si are highly depleted in dust; see Li & Greenberg 1997 and 
references therein); [X/HL t be the amount of X relative H locked up in dust ([C/HL t = 
[C/H] -[C/H] gas « 251 ppm, [N/H] dust « 24.2 ppm, [0/H] dust « 235 ppm, [Mg/H] 

dust ~ 34.5 ppm, 

[Fe/H] dust ~ 34.4 ppm, [Si/H] dust Ri 28.1 ppm). Assuming a stoichiometric composition of MgFeSiC>4 
for interstellar silicates, the total mass of silicate dust per H atom is m s n « [Fe/H] dust /j,f c + 
[Mg/H] dust fi Mg + [Si/H] dust usi + [0/H] sil no ss 5.61 x 10~ 3 /in where fix is the atomic weight of X 
in unit of m « 1.66 x 10" 24 g, and [0/H] sil « 4 ([Fe/H] dust + [Mg/H] dust + [Si/H] dust )/3 « 129 ppm 
is the amount of O in silicate dust per H atom (i.e., we assign 4 O atoms for the average of the 
Fe, Mg, and Si abundances). The carbonaceous dust component is dominated by C, with lit- 
tle H, N, and O (we assume H/C=0.5, O/C=0.1). The total mass of carbon dust per H atom 
is m carb « [C/H] dust/ x c + [N/H] dustm + 0.5 [C/H] dust m + 0.1 [C/H] du8t fi « 3.88 x 10~ 3 m . 
The C, O, and N atoms left over after accounting for the silicate and carbon dust components 
are assumed to condense in icy grains in the form of H2O, NH3, CO, CO2, CH3OH and CH4 
(following Greenberg [1998], we assume CO:CO 2 :CH 3 OH:CH4:H 2 CO=10:4:3:l:l). The total mass 
of icy grains per H atom is m^ cc ~ m^ e + mf ce + m™ ter , where the mass of C-containing ice 
mg e w [C/H] gas ^ c + [C/H] gas (22no + 18/x H )/19 w 2.87 x lO" 3 ^; the mass of NH 3 ice mf ce w 
[N/H] gas ( MN + 3^ H ) ~ 4.10 x 10~Vh; the mass of water ice m™ tcr w [0/H] watcr (// Q + 2 m ) w 
4.12 x 10- 3 ^ H ; [0/H] watcr « [O/H] - [0/H] sil - 0.1 [C/H] dust - 22 [C/H] gas /19 « 229 ppm is the 
amount of O locked up in H2O ice (we assume H2O contains all the remaining available O). 



7 Some H will be present, mostly in combination with O, C, and N. 
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Therefore, as a first approximation, we may assume a mixing ratio of m car b/w- s ii ~ 0.7 and 
^-carb) ~ 0-8 for cold regions (for hot regions where ices sublimate the dust can be 
simply modelled as porous aggregates of silicate and carbon particles with m car h/m s n ~ 0.7). 
This does not deviate much from the in situ measurements of cometary dust (m car b/m s ii ~ 0.5, 
mi cc /[m s ii + m car b] ~ 1.0; see Greenberg & Li 1999 and references therein) which is often suggested 
to be porous aggregates of unaltered interstellar dust (Greenberg 1998; Greenberg & Li 1999). 

B. Dust Morphology 

Dust formed through coagulation of many small subgrains has a porous structure. Let v s n 
and f car b be the total volumes taken up by the silicate and carbonaceous dust components, re- 
spectively; P be the porosity - the fractional volume of vacuum; p s ii (~ 3.5 gem -3 ) and p car b (~ 
1.8 gem -3 ) be the mass densities of silicate and carbonaceous materials, respectively. The mass 
density for a fluffy aggregate of silicate and carbonaceous subgrains with a porosity of P is 
(P) = (l-^)(p s ii«sii + /'carbWcarb)/(vsii + Vcarb)- In regions colder than « 110 - 120 K, pre- 
existing or recondensed ices around the silicate and carbonaceous dust cores would fill all or part 
of the vacuum in the aggregate. Let v- lcc be the total volume taken up by the ice component; 
Pice (~ 1-2 g cm -3 ) be the ice mass density. The porosity for the fluffy aggregate of ice-coated silicate 
and carbonaceous subgrains would be reduced to P' = max {0, 1 — (1 — P) [1 + v- lcc / {v s \\ + f ca rb)]}- 

Its mean mass density is (//) = (1 - P') (p s i\V si i + Pcarb^carb + Pice Vice) / (v si \ + ^carb + Vice)- The 

dust "size" is defined as the radius of the sphere which encapsulates the whole fluffy aggregate. 
The mean mass densities defined here are independent on aggregate size. 

REFERENCES 

Augereau, J.C., Lagrange, A.M., Mouillet, D., Papaloizou, J.C.B., & Gorod, P.A. 1999, A&A, 348, 
557 

Backman, D.E., & Paresce, F. 1993, in Protostars and Planets III, ed. E.H. Levy & J.I. Lunine 
(Tucson: Univ. Arizona Press), 1253 

Beckwith, S.V.W., Henning, Th., & Nakagawa, Y. 2000, in Protostars and Planets IV, ed. V. 
Mannings, A. P. Boss, & S.S. Russell (Tucson: Univ. Arizona Press), 533 

Bohren, C.F., & Huffman, D.R. 1983, Absorption and Scattering of Light by Small Particles (New 
York: Wiley) 

Burns, J. A., Lamy, P.L., & Soter, S. 1979, Icarus, 40, 1 
Draine, B.T., & Lee, H.M. 1984, ApJ, 285, 89 
Gail, H.-P., 2001, A&A, 378, 192 



16 



Greenberg, J.M. 1998, A&A, 330, 375 

Greenberg, J.M., k Li, A. 1999, Space Sci. Rev., 90, 149 

Jayawarelhana, R., Fisher, R.S., Hartmann, L., Telesco, CM., Piha, R.K., k Fazio, G. 1998, ApJ, 
503, L79 

Jura, M. 1991, ApJ, 383, L79 

Jura, M., Zuckerman, B., Becklin, E.E., k Smith, R.C. 1993, ApJ, 418, L37 

Jura, M., Ghez, A.M., White, R.J., McCarthy, D.W., Smith, R.C, k Martin, P.G. 1995, ApJ, 445, 
451 

Jura, M., Malkan, M., White, R., Telesco, CM., Piha, R.K., k Fisher, R.S. 1998, ApJ, 505, 897 
Kenyon, S.J., Wood, K., Whitney, B.A., k Wolff, M.J. 1999, ApJ, 524, L119 
Klahr, H.H., k Lin, D.N.C. 2000, ApJ, 554, 1095 

Koerner, D.W., Ressler, M.E., Werner, M.W., k Backman, D.E. 1998, ApJ, 503, L83 

Kurucz, R.L. 1979, ApJS, 40, 1 

Li, A., k Draine, B.T. 2001a, ApJ, 550, L213 

Li, A., k Draine, B.T. 2001b, ApJ, 554, 778 

Li, A., k Greenberg, J.M. 1997, A&A, 323, 566 

Li, A., k Greenberg, J.M. 1998, A&A, 331, 291 

Li, A., k Lunine, J.I. 2003, ApJ, 594, 987 

Schneider, C, et al. 1999, ApJ, 513, L127 

Shu, F.H., Shang, FL, k Lee, T. 1996, Science, 271, 1545 

Sitko, M.L., Lynch, D.K., k Russell, R.W. 2000, AJ, 120, 2609 

Sofia, U.J., k Meyer, D.M. 2001, ApJ, 554, L221 

Spangler, C, Sargent, A.I., Silverstone, M.D., Becklin, E.E., k Zuckerman, B. 2001, ApJ, 555, 932 
Telesco, CM., et al. 2000, ApJ, 530, 329 

Weidenschilling, S.J., k Cuzzi, J.N. 1993, in Protostars and Planets III, ed. E.H. Levy k J.I. Lunine 
(Tucson: Univ. Arizona Press), 1031 

Wooden, D.H. 2002, Earth, Moon, k Planets, 89, 247 



17 



Wyatt, M.C., Dermott, S.F., Telesco, CM., Fisher, R.S., Grogan, K., Holmes, E.K., & Piha, R.K. 
1999, ApJ, 527, 918 

Zuckerman, B., & Becklin, E.E. 1993, ApJ, 414, 793 



This preprint was prepared with the AAS IATjrjX macros v5.0. 



18 



1 1 1 1 1 1 I I: 

HR 4796A 

cold -coagulation 
dust model ~ z - 




j i i i i i i i 



30 100 

A (/xm) 



1000 



Fig. 1 . — Comparison of the observed spectral energy distribution of the HR 4796A dust disk to theoretical IR emis- 
sion spectra calculated from the "cold-coagulation" model (porous aggregates of amorphous silicate dust and carbona- 
ceous dust plus ices in cold regions). The dust in this ring-like disk is taken to have a Gaussian distribution function 
a{r) = cr p exp[-4 In 2{(r -r )/A} 2 ] with its peak at r p = 70 AU and a FWHM A = 15 AU (r € [0.15 AU, 250 AU]) 
as inferred from the scattered light images (Schneider et al. 1999; Kenyon et al. 1999). We assume a power-law dust 
size distribution dn(a)/da oc a~ a (a £ [1 /im, 1cm]). Upper panel (a): dotted line (model no. 1) - best-fit P = 0.95 
(P' « 0.87; see §B) model (a w 2.8, m d « 2.67 x 10 27 g, a p « 2.97 x 10 4 cm" 2 , yf/N « 2.85 where N = 18 [16 
data points plus 2 free parameters P and a; see text]); solid line (model no. 2) - best-fit P = 0.90 (P' ~ 0.73) model 
(a w 2.9, m d ~ 4.03 x 10 27 g, <r p w 4.90 x 10 4 cm" 2 , x 2 /N ~ 1-81). Lower panel (b): dotted line (model no. 3) - 
best-fit P = 0.80 (P' w 0.46) model (a w 3.0, m d w 5.73 x 10 27 g, a p w 7.55 x 10 4 cm" 2 , x 1^ ~ 6.08); solid line 
(model no. 4) - best-fit P = 0.60 (P' w 0) model (a w 3.2, m d w 4.34 x 10 27 g, cr p w 1.31 x 10 5 cm -2 , x /N w 7.60). 
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Fig. 2. — Same as Figure 1 but for the "hot-nebula" model (porous aggregates of crystalline silicate dust plus 
ices in cold regions). Upper panel (a): dotted line (model no. 5) - best-fit P = 0.95 (P' w 0.90) model (a « 2.6, 
m d « 6.02 x f0 27 g, <7p « 8.54 x 10 3 cm" 2 , % 2 /iV « 5.29); solid line (model no. 6) - best-fit P = 0.90 (P' w 0.80) 
model (a w 2.8, m d w 7.48 x 10 27 g, cr p w 2.55 x f0 4 cm~ 2 , \ 2 / N ~ 4 -48). Lower panel (b): dotted line (model 
no. 7) - best-fit P = 0.80 (P' w 0.60) model (a w 2.9, m d w 1.16 x 10 28 g, cr p w 4.30 x 10 4 cm" 2 , xV^ « 5.24); 
solid line (model no. 8) - best-fit P = 0.60 (P' w 0.20) model (a » 3.1, m d w 9.25 x 10 27 g, cr p w 8.04 x 10 4 cm" 2 , 
X 2 /N w 9.94). 
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Fig. 3. — Comparison of the HR4796A observational SED to theoretical IR emission spectra calculated from (a) the 
"cold-coagulation" model and (b) the "hot-nebula" model. The dust spatial distribution is assumed to be a power-law 
a(r) = er p (V/Vp) - * 3 for either 0.15 < r < 250 AU (dotted lines) or 40 < r < 130 AU (solid lines). We assume a power 
law dust size distribution dn(a)/da oc a~ a (a G [1 /im, 1cm]). Upper panel (a): dotted line (model no. 9) - best-fit 
P = 0.90 (P' w 0.73) "cold-coagulation" dust model for r G [0.15, 250 AU] (0 w -2.3, a w 3.2, m d w 8.43 x 10 2r g, 
cr p w 4.17 x 10 3 cm -2 , xV^ ~ 2.32); solid line (no. 10) - best-fit P = 0.90 (P' w 0.73) "cold-coagulation" dust 
model for r G [40, 130 AU] (/3 w -2.5, a w 3.0, m d w 5.81 x 10 27 g, cr p w 9.75 x 10 3 cm" 2 , x 2 /^ w 3.28). Lower 
panel (b): dotted line (model no. 11) - best-fit P = 0.90 (P' w 0.80) "hot-nebula" dust model for r G [0.15, 250 AU] 
(/3 w -1.1, a w 2.9, m d w 2.55 x 10 28 g, <r p w 2.58 x 10 3 cm -2 , x 2 /N w 9.71); solid line (model no. 12) - best- 
fit P = 0.90 (P' w 0.80) "hot-nebula" dust model for r G [40, 130 AU] (f3 w 5.2, a w 2.6, m d w 6.29 x 10 27 g, 
cr p w 5.55 x 10 2 cm" 2 , X 2 /N « 2.02). 
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Fig. 4. — Comparison of the HR4796A observational SED to theoretical IR emission spectra calculated from the 
"cold-coagulation" model (solid line [no. 14]; P = 0.90, P' w 0.73, a w 2.9, m d w 5.41 x 10 27 g, cr p » 2.34 x 10 4 cm" 2 , 
X 2 /N » 3.60) and the "hot-nebula" model (dotted line [no. 13]; P = 0.90, P' w 0.80, a w 2.8, m d » 1.00 x 10 28 g, 
(T P w 1.22 x 10 4 cm" 2 , x 2 / N « 4.89). The dust spatial distribution is taken to be that of Telesco et al. (2000). 
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Fig. 5. — IR emission from dust in the ring regions of n n = 0.15 AU < r < r out = 250 AU (solid line; same as the 
solid line in Figure la), (r p — A) < r < (r p + A) (dotted line), (r p — A) < r < r p (dashed line), and r p < r < (r p + A) 
(long-dashed line) calculated from the best-fitting "cold-coagulation" dust model (no. 2; P — 0.90, P' w 0.73, a w 2.9) 
with a Gaussian spatial distribution (r p = 70 AU, A = 15 AU). The (r p — A) < r < (r p + A) spectrum (dotted line) 
is indistinguishable from the 0.15 AU < r < 250 AU spectrum (solid line). 



23 



I 



10- 9 

10 -io 

10- 11 
io- 12 b 



e i — r~n 



7 10- 13 = 



CO 
Clfl 

S-H 

CD 



10- 



10 



10 _ 



< IO" 11 



IO" 12 

io- 13 



— I 1 1 — I — I — I I 1=1 

HR 4796A 

@\ cold-coagulation 
dust model 



- / 




r p =70 AU,\A=15 AU 
a(r) = CT p exp[-41n2.[(r-r p )/A] 2 ] 
= 0.90, P'=0.73, (T p =4:90xl0 4 cm- 2 

ljum < a < 1 cm 

ljum < a < 10/xm 

— < a < 100/zm 



10 



30 100 

A (/xm) 



300 



1000 



Fig. 6. — Upper panel (a): models with smaller a m i n (= 0.1 /im; solid line [no. 15]; a « 2.8, « 5.49 x 10 27 g, 
Op w 1.93 x 10 6 cm- 2 , x 2 /^ « 3.44) and larger o min (10 /an; dotted line [no. 16]; a « 3.1, m d » 2.98 x 10 27 g, 
a p ~ 1.35 x 10 3 cm -2 , x 2 /N ~ 5.88). Other model parameters are the same as those of the best-fit "cold-coagulation" 
model (no. 2; P — 0.90, P' ~ 0.73). Lower panel (b): model spectra calculated from the best-fit "cold-coagulation" 
model (no. 2; P = 0.90, P' w 0.73, a w 3.1) with 1 /im < a < 10 ^m (dotted line), 1 /im < a < 100 /xm (dashed line), 
1 /xm < a < 1 cm (solid line; same as the solid line in Figure la). 
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Fig. 7. — Upper panel (a): IR emission per gram dust mass (received at the Earth) calculated from the P = 0.90 
(P' ~ 0.73) "cold-coagulation" dust at r = 70 AU with a single size of a = 1/xm (solid line), 10pm (dotted line), 
100 /im (dashed line), 1000 fim (long-dashed line), and 1cm (dot-dashed line). Lower panel (b): IR emission per 
gram dust mass calculated from the P = 0.90 (P' fa 0.73) "cold-coagulation" dust at r — 55 AU (dotted line), 70 AU 
(solid line), 85 AU (dashed line), 130 AU (long-dashed line), and 250 AU (dot-dashed line) with a distribution of sizes 
(dn/da oc a~ a , a — 2.9, a £ [1 fim, 1 cm]). 
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Fig. 8. — Upper limits on the "zodiacal" dust component (a) at r = 9AU (ma ~ 9.5 x 10 -6 me) and (b) at 
r = 4.5 AU (rrid ~ 1.8 x 10 -6 me): dashed lines - model spectra calculated from the "zodiacal dust" which is taken to 
be porous crystalline silicate dust with P = 0.97 and a — 450 /im (see Augereau et al. 1999) of (a) rrid ~ 9.5 x 10 -5 mj 
at r — 9 AU and (b) nid w 1.8 x 10 -5 mj at r — 4.5 AU, assuming the emission at A < 12 /im is entirely produced by 
the "zodiacal" dust; dotted lines - the best-fit "cold-coagulation" model spectra (no. 2; P = 0.90, P' ~ 0.73, a w 2.9; 
same as the solid line in Figure la); solid lines - the sum of 10% of the "zodiacal" dust spectrum and the best-fit 
"cold-coagulation" model spectrum. 
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Fig. 9. — Upper panel (a): the ratios of the radiative repulsion to the gravitational attraction (/3rp) for the 
best-fitting "cold-coagulation" dust (without ice P — 0.90 [solid line] or with ice P' — 0.73 [dotted line]) and the 
"hot-nebula" dust (without ice P = 0.90 [dashed line] or with ice P' = 0.80 [dot-dashed line]). The long-dashed 
horizontal line plots /3rp = 1. Lower panel (b): the orbit decay timescales tpr due to the Poynting-Robertson drag 
for the best-fitting "cold-coagulation" dust (without ice P = 0.90 [solid line] or with ice P' = 0.73 [dotted line]) 
and the "hot-nebula" dust (without ice P — 0.90 [dashed line] or with ice P' = 0.80 [dot-dashed line]) at a radial 
distance of r = 70 AU from the central star (note tpr oc r 2 ). The long-dashed horizontal line plots the HR4796A age 
(«8Myr). 



